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Abstract. We explore the consequences of the excluded volume interaction of hard 
spheres at high densities and present a theory for excited granular materials. We 
first demonstrate that, in the presence of gravity, the granular density crosses over 
from Boltzmann to Fermi statistics, when temperature is decreased in the weak 
excitation limit. Comparisons of numerical simulations with our predictions con- 
cerning the scaling behavior of temperature with agitation frequency, gravity and 
particle-diameter show satisfying agreement. Next, within the framework of the En- 
skog theory of hard spheres, we interpret this crossover as a "condensation" of hard 
spheres from the dilute gas-state to a high density solid-like state. In the high den- 
sity, low temperature limit Enskog theory fails because it predicts densities larger 
than the closed packed density below a certain temperature. We show how to ex- 
tend the range of applicability of the Enskoq theory to arbitrarily low temperatures 
by constructing a physical solution: all particles that are situated in regions with 
densities larger than a certain maximum density are assumed to be "condensed" . 



1 Introduction 

This paper is a review of a recently proposed theory of granular dynam- 
ics [0J|] , which is based on the simple recognition that granular materials are 
basically a collection of hard spheres that interact with each other via a hard 
core potential H . For this reason, many of the properties of excited granular 
materials may be understood from the atomistic view of molecular gases, in 
particular from the point of view of kinetic theory Q. There are, however, 
several distinctions between molecular gases and granular materials: First, 
the mean free path of the grains can be rather small - less than the particle 
diameter - even if the system is strongly forced. Second, granular material 
is made of macroscopic particles with finite diameter, and thus the mate- 
rial cannot be compressed indefinitely. When the mean free path vanishes, 
the density is the maximum, closed packed density. Third, gravity plays an 
important role in the collective response of granular materials to external 
stimuli because of the ordering of the particles according to their potential 
energy in the gravitational field. 

For example, one of the notable characteristics of excited granular mate- 
rials in a confined system under gravity is the appearance of a thin boundary 
layer near the surface that separates a fluidized region from a solid region. 
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This effect is also seen in shearing experiments [|7j , avalanches || , and grains 
subjected to weak excitations P,^p[. In this limit, those grains in a solid re- 
gion are effectively frozen, and thus do not participate in dynamical, diffusive 
processes. Hence, the conventional Boltzmann statistics, which is applicable 
in the limit of strong excitation and rapid flow where all the particles are 
dynamically active, certainly needs modification. Our first aim is to show 
that, in the weak excitation limit, the statistics where such a boundary layer 
appears, is analogous to the Fermi statistics. We will in fact demonstrate 
that the density profile of the grains is qualitatively well given by a Fermi 
function, from which we define the global temperature T and develop a ther- 
modynamic theory of configurational statistics for excited granular materials. 
We present an explicit formula to relate the temperature T to the external 
control parameters such as the frequency u> and amplitude A of the excitation 
vibration, the diameter D of the grains, and the gravitational acceleration g. 
Next, we examine the microscopic basis of the crossover from Boltzmann to 
Fermi statistics based on Enskog theory for hard spheres jnj, and demon- 
strate how the crossover proceeds as grains condense from the bottom toward 
the surface. 



2 Configurational Statistics and Maximum Entropy 
Principle: Justification of a Thermodynamic 
Approach 

It is quite well known that variational methods are not useful in determining 
the properties of nonequilibrium systems. In systems that show e.g. periodic 
patterns, it is evident that one cannot perform variations about a steady 
state situation. Since the system being studied here is a dense, dissipative, 
nonequilibrium granular system, we find it necessary to make some comments 
on this point. If the mean free path of the grains is much less than a particle 
diameter, each particle may be considered to be effectively confined in a cage 
as also assumed for the free volume theory of dense liquids [Q. In such 
a case, the basic granular state is not a gas, but a solid or a crystal |J. 
Thus an effective thermodynamic theory based on the free energy argument 
may be more appropriate than the kinetic theory. Our argument is that the 
dense state can be assumed as a "steady state" for which we compute the 
"configurational statistics" by means of the usual variational method as the 
most probably state. 

To be more specific, consider the excitation of disordered granular mate- 
rials confined in a box with vibrations of the bottom plate. The vibrations 
will inject energy into the system which cause the "ground state" to become 
unstable. A new, excited state will emerge with an expanded volume. The 
time average of such configurations which have undergone structural distor- 
tions, may be deviating weakly from the ground state so that the use of an 
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effective thermodynamic theory based on the variational principle could be 
justified. 

Such a thermodynamic approach may be further justified by the following 
two recent experiments in both the weakly and the strongly excited regimes: 

• Weakly or moderately excited regime: Clement and Rajchenbach(CR) || 
have performed an experiment with the vibrational strength, i.e. the di- 
mensionless vibration acceleration r = Auj 2 /g, for a two dimensional 
vibrating bed, using inclined side walls to suppress convection. CR have 
found that the ensemble-averaged density profile as a function of height 
from the bottom layer obeys a universal function that is independent of 
the phase of oscillations of the vibrating plate. Namely, it is independent 
of the kinetics imposed on the system. One conceptually important point 
here is that the reference point of the density profile is not the bottom 
plate, but the bottom layer, which of course is fluidized. 

• Strongly excited regime: Warr and Hansen (WH) || have performed an 
experiment on highly agitated, vertically vibrating beds of r ps 30 — 50 
using steel balls with a large coefficient of restitution. They have found 
that the collective behavior of this vibrated granular medium in a sta- 
tionary nonequilibrium state exhibits strong similarities to those of an 
atomistic fluid in thermal equilibrium at the corresponding particle pack- 
ing fraction, in particular, concerning the two-point correlation function 

§■ 

The results of both experiments indicate that for both moderately and 
strongly excited systems, a one-to-one correspondence seems to exist between 
the configurational statistics of the nonequilibrium stationary state and the 
equilibrium thermal state. In fact, this is not so surprising considering the fact 
that upon vibration, the granular materials expand and consequently the vol- 
ume of the system increases. In turn, this increase corresponds to a rise in 
the potential energy after the configurational average is appropriately taken. 
Then the problem reduces to the packing problem, and the temperature-like 
variable, T, may be associated to the vibrating bed. The existence of distinc- 
tive configurational statistics in the density profile of CR (and also in WH in 
a special case) appears to be fairly convincing evidence that kinetic aspects 
of the excited granular materials may be separated out from the statistical 
configurations. These observations are the basis of the thermodynamic theory 
proposed in [0. 

3 Fermi Statistics of Weakly Excited Granular 
Materials 

We first view the system of granular particles as a mixture of holes and 
particles as in the lattice gas or the diffusing void model |l2|, which is the 
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simplest version of the free volume theory Jl3| . We now assign virtual lattice 
points by dividing the vibrating bed of width L and the height /j,D, with D 
typically the diameter of a grain and \x the number of layers, into cells of size 
D x D. Each row, i, is then associated with the potential energy ti — mgzi 
with Zi = (i — 1/2)1) and m the mass of the grain. Note that the degeneracy, 
fl, of the each row is simply the number of available cells, i.e: Q = L/D, For a 
weakly excited system with r ~ 1, the most probable configuration should be 
determined by the state that maximizes the entropy in the micro-canonical 
ensemble approach. 

Taking into account the excluded volume effect which plays a similar role 
for dense granular systems as the Pauli principle in Fermi statistics, we derive 
the entropy S, defined as the total number of ways of distributing N particles 
into the system. Standard counting argument E3] yields, 



where iV, is the number of particles occupying the i—th row. Since gravity 
orders the grains with respect to their potential energy, a grain can be seen 
as a "spinless Fermion" , where the height z plays the role of the momen- 
tum variable (if one makes the connection to the electron gas). Maximizing 
S with the constraints, Y,i N * = N and Y,i N i e i = ( U (T(r))), the mean 
steady state system energy, we find that the density profile, 4>(z), which is 
the average number of occupied cells at a given energy level, is given by the 



where f3 = mgD/T, the height is z — zjD and analogous to the Fermi 
energy, fx, measured in units of D, is the number of layers in the low tem- 
perature limit, i.e. in the system at rest. Note that both \i and T enter as 
Langrange multipliers introduced by the two constraints, i.e, the conserved 
number of particles and the mean energy. The global temperature T defined 
here is similar to the compactivity introduced by Edwards and his collabo- 
rators in their thermodynamic theory of grains |l6| , but is different from the 
kinetic temperature defined through the kinetic energy JT7| . We point out 
that the Fermi statistics is essentially the macroscopic manifestation of the 
classical excluded volume effect and the anisotropy which causes the ordering 
of potential energy by gravity. In the spirit of the proposed analogy, the top 
surface of the granules at rest plays the role of the Fermi surface, and the 
grains in the thin boundary layer near the surface play the role of the excited 
electrons of a Fermi gas in metals. 

For strongly excited systems, the exclusion principle does not apply. The 
Fermi analogy is valid when the zero temperature Fermi energy satisfies [i 3> 
ni, where n; is the number of fluidized layers. Now, the energy £7j, injected 
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Fermi distribution 



<j>{z) = N./Q = 1/[1 + cxp(/3(z - n))\, 



(2) 
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into the system is of the order of mA 2 Lo 2 jl and the potential energy E p , 
needed to fLuidize the particles in the top ni layers is of the order of mgniD. 
Equating these two energies we find a necessary condition for the fluidization 
of the top ni layers, namely: ni ~ rA/(2D). Hence, we expect the Fermi 
statistics to be valid for P ~ 1, if fi 3> A/D, 

4 Relation between T and J 1 

We now relate the temperature T to the external control parameters such as 
r. First, a thermal expansion. We determine the energy per column u{T) = 
J 4>(z)mzgdz, from which we can determine the shift in the center of mass 
per particle; h{T) = u(T)/mg, which is given by: 



with h(0) = Dii 2 /2. Second, a kinetic expansion. We make an observation 
that for a weakly excited granular system, most excitations occur near the 
Fermi surface, and thus the volume expansion may be effectively well repre- 
sented by the maximum height, Ho(r) of a single particle bouncing on the 
vibrating plate assuming that the Fermi surface is in contact with the vibrat- 
ing plate. The kinetic expansion, H (r), is then determined by the maximum 
of A(t) in the following equation that describes the trajectory of a single ball 
on a vibrating plane with intensity r = oJ 2 /g with A the amplitude, u> the 
frequency of the vibrating plate, and g the gravitational constant: 



in units of g = uj = 1, where t = sin _1 (l// 1 ). Note that since A is the relative 
distance between the ball and the plate, it cannot be negative. More precisely, 
the particle is launched from the plate at t = t by inertia, and then makes 
a free flight motion until it strikes the plate. It then stays on the plate until 
t = t + 2n, when it is launched again and repeats the same motion. Hence, A 
is a periodic function with period 2-k. Since the Fermi distribution near T = 
can be approximated by a piecewise linear function and Hq (r) is thought to 
be the edge of the function, we expect H (r) w Ah/2 = (h(T) - h(0))/2. By 
equating the thermal expansion, (||), to the kinetic expansion, Ho(r)g/uj 2 , in 
physical units, we now complete our thermodynamic formulation by present- 
ing the explicit relation between T and _T Q: 



In MD simulation, one may measure the maximum height, h a (r), of a single 
ball on a vibrating plate and replace gH {r)/uj 2 with 2h (r)/a with a an 




(3) 



A(t) = r(sm(t ) - sin(t)) + rcos{t )(t - to) - -(t - t ) 2 (4) 




adjustable parameter Q. 
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Now we compare our theoretical prediction with an experimental result of 
Clement and Rajchenbach Figure 1 shows the fitting of the experimental 
density profile for F — 4 of CR by the scaled Fermi distribution, <f>(z) = 
p(z)/p c , with p c the closed packed density. For hexagonal packing, p c w 0.92. 
The fitting value of T/mg is 2.0 mm, while Eq.(|) yields T/mg m 2.6 mm. The 
agreement between the two is fairly good in spite of such a simple calculation. 
This expression also agrees with the simulation result p]Jl5| ■ 
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Fig. 1. Density <f>(z) as a function of the height z. The symbols are the data by 
Clement and Rajchenbach and the dotted line is the Fermi distribution function, 
Eq.(2). 



Note that the detailed expression of H (T) depends on the manner by 
which the grains are excited and we expect that our main scaling prediction 
of Eq.(||), namely T cx g z l 2 D 1 l 2 /ui, will hold even for systems driven not by 
sinusoidal waves. Further, T has a gap at T = because the time between 
the launching and landing of the ball is always finite for F > 1. Next, it is 
well known that the specific heat per particle, C v = du/dT, can be written 
as the fluctuations in the energy, namely ((Au) 2 ) — ((u(z) — u) 2 ) — T 2 C V 
| fl8| . Hence, our theory makes a nontrivial prediction for the fluctuations in 
the center of mass: 

«^> = <«r)^> = «^4(^) 3 ^ ,6) 

Mo 3 \mgDJ \i% 

while the center of mass is given by: 

Az(T) = z(T) - z(0) = ^ ( -I—) 2 (7) 

Note that the total expansion, Ah(T) = p Az and its fluctuations ((Ah) 2 ) /D 2 
= (no(Az) 2 / D 2 ) are only a function of the dimensionless Fermi temperature 
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Tf = T/mgD as expected. Furthermore, note that (6) is an indirect confir- 
mation that the specific heat is linear in T as it is for the non-interacting 
Fermi gas. 

5 Test of Fermi Statistics by Molecular Dynamics 
Simulations 

In this section, we examine the configurational statistics of granular materials 
in a vibrating bed, in particular the density profile, and its fluctuations by 
MD simulations and compare the results with the predictions made in the 
previous chapter. The MD code was provided by the authors of and the 
details of the code can be found in the literature. Simulations were carried 
out in two dimensional boxes with vertical walls for N particles each with 
a diameter D = 0.2 cm and a mass m — 47r(Z?/2) 3 /3 with the degeneracy 
Q using (JV, Q) = (100,4), (200,4) and (200,8) with a sine wave vibration. 
The dimensionless Fermi energy /j, = N/Q is the system height at rest. For 
all cases, the inequality /i ^> A/D was satisfied. In Fig. 2 the temperature, 
obtained by the best fit of the density profile to the Fermi function (dots), is 
plotted against r and the values predicted by Eq.(5). Note the fairly good 
agreement between theory and simulations. 

We also studied the temperature scaling against the frequency, gravity 
and diameter to further check the validity of Eq. (5). The scaling laws as 
determined by the simulations are: 





i oj-" 11 


T', 




T p 





with mi ps 1.16, TO2 ~ 0.48, and ps 0.53. These values are close to the 
predicted ones by Eq. (5), i.e m\ = 1, mi = 0.5, = 0.5. For detailed 
comparisons, see the original paper |l5| . 

Finally, we have also checked the scaling of the center of mass, and the 
fluctuations against T 2 and T 3 . Since the density profiles are well fitted by 
the Fermi function, we anticipate that the center of mass and its fluctuations 
obey the scaling as shown in Figs. 3 and 4. Note that the increase in the 
center of mass is second order in temperature T, which is contrary to the 
mean field prediction of the linear increase of volume in the compactivity X 

6 Condensation of Hard Spheres Under Gravity 

Our next goal is to examine the microscopic basis of Fermi statistics based 
on the kinetic theory, in particular the Enskog equation for hard spheres of 
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Fig. 2. Comparison between the measured temperatures by MD(dots) and the 
predicted ones. The lower line was obtained by Eq.(5), while the upper one was 
obtained by replacing gH /uj 2 by 2h (r)/a with a w 12.8. h (r) is the maximum 
jump height of a single ball on a vibrating plated obtained by MD. 

mass m and diameter D, to explore whether or not the kinetic theory can 
describe the cross over from Boltzmann to Fermi statistics and if so, under 
what conditions it occurs. 

Our particularly interesting discovery Q is that the prediction of the En- 
skog equation is only valid when (3/i < /i , where /i is the dimensionless initial 
layer thickness of the granules (or the Fermi energy), (3 — mgD/T with T 
the temperature, and the critical value, \i 0l is determined to be \i Q — 21.756 
in two dimensions (2D) and /i = 15.299 in three dimensions (3D). (For max- 
imum packing, /i Q = 151.36 in 3D and [i a = 144.6155 in 2D.) When this 
inequality is violated, the Enskog equation does not conserve the particles, 
and we present a scenario of resolving this puzzle based on physical intu- 
ition, namely that the missing particles condense from the bottom toward 
the surface In this way, the hard sphere Enskog gas appears to contain 
the essence of Fermi statistics. Wc briefly describe how the density profile, <fi 
as a function of dimensionless variable can be obtained from the Enskog 
equation. For details, see Ref. [2]. Note that (j> introduced here is not the 
volume fraction v but is scaled so that <fi = 1 at closed packing. 

In a free volume theory, particles are confined in a cage. Hence, if we 
use a simple cubic lattice as a basic lattice, the close packed volume fraction 
p c = N/V = N/(D 2 N) = 1/D 2 . If we define the dimensionless density 
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Fig. 4. Fluctuations in the center of mass as a function of temperature. Symbols 
used are the same as in Fig. 3. 
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<j>(z) = G{z)/p c = D 2 G{z) or 4>((,/3) = D 2 G{z) with C = z/D, we then 
obtain the following exact dimensionless equation of motion for (/>(£,/?): 



where = mgD/T and 



</ " ( ° ' .%(()= 0(C)/ C (C) (8) 



1 f 2 * .„ r I I 



^C(C) = 2 / d0cos0[x(C- 2Cos0MC-cos0)-x(C+2Cos0)«KC + cos0)] 

(9) 

For 3D, the corresponding equation for the density <p(() = D 3 G(z) is given 
by Eq.(8) with: 



I C (Q =7r f dOsinO cos 



x(C~cos 6»/2)0(C-cos 6»)-x(C+cos 6»/2)0((+cos ( 



(10) 

Several forms for the equilibrium correlation function \ have been proposed, 
but we use the following widely used forms: For 2D, we use the form proposed 
by Ree and Hoover |^0): x{4>) = (1 - "10 + "2</> 2 )/((l — acf) 2 with a = 
0.489351 • tt/2 « 0.76867, a x = 0.196703 • tt/2 w 0.30898, a 2 = 0.006519* 
7r 2 /4 w 0.0168084, while for 3D, we use the form suggested by Carnahan and 
Starling |§: X {4>) = (1 - ^/12)/(1 - ^/6) 3 

Since the total number of particles, iV, remains fixed, the following nor- 
malization condition should be satisfied for both 2D and 3D. 

dC0(C;/3) = M (ii) 

where fi = N / Q x (or /i = N / ft x f} y in 3D) is the Fermi energy and fl x , £2 y 
are the degeneracies along the x and y axes. We now perform the gradient 
expansion of (9) and (10) and retain only the terms to first order in d\j 'd£. 
The 3D solutions for the first order differential equation can be obtained 
easily, and are given by || : 

- P(C - p) = ln4> - 1/(1 - acf) 2 + 2/(1 - a(t>f (12) 
pji = ln(4> ) - 1/(1 - a^) 2 + 2/(1 - act> f (13) 

= - -. — + jz —-^ (14) 

l-a(p (l-acpo) 6 

where a = it/6. For given values of (3 and /i, cj> = 0(£ = 0) will be determined 
by Eq.(14) with the condition that <f> cannot in anyway greater than the 
closed packed density. As we will see, we need care to proceed. First, since 
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the right hand side of (14) is an monotonically increasing function of <j> and 
<f> cannot be greater than the closed packed density, which is 1 in our units, 
(3p must have an upper bound p ; namely, p = 15.299 in 3D, (p D = 21.756 
in 2D), which is the value obtained by setting <p = 1 in the right hand side 
of (14). Note that this upper limit depends on the underlying basic lattice 
structure. (For a hard sphere gas in a continuum space, at the closed packed 
density, r) = ^(f ) 3 = f <j> w 0.74 in 3D and r\ = \ = tt/2\/3 w 0.907 in 
2D hexagonal packing, in which case, the upper limit p a = 151.36 in 3D and 
p Q = 144.6155 in 2D.) Considering the fact that both the temperature T and 
the Fermi energy p, are arbitrary control parameters, the existence of such 
bounds is a puzzle: if ftp is less than p a , then the density profile given by 
Eq.(12) is well defined, but if ftp is greater than p„, then <j> must be one at 
the bottom, and the particle conservation appears to break down, namely 



The central question is: where does the rest of the particles go? In or- 
der to gain some insight into this question, consider first the case of point 
particles under gravity, in which case the density profile is given by: p(() = 
p(0)exp(—mg(/T). If we put more particles into the system, we simply need 
to increase p(0) because the point particles can be compressed indefinitely, 
and the profile simply shifts to the right. We now replace these point par- 
ticles with hard spheres, which cannot be compressed indefinitely and thus 
the maximum density at any point is the closed packed density. Suppose we 
have a system of hard spheres at a certain temperature T, where the density 
is closed packed at the bottom layer and smoothly decreases toward the sur- 
face. At this point, if we add more hard spheres, say the amount of one layer, 
to the system, how does the density profile modify? Since the hard spheres 
cannot be compressed, our intuition tells us that after the system reaches the 
equilibrium, the density of the first (bottom) and second layer becomes closed 
packed, forming a rectangle, which is then followed by the original Enskog 
profile. If we keep adding more particles, we obtain the density profile that 
is the combination of the rectangle (we term this the Fermi rectangle) be- 
ginning at the bottom and the smooth original Enskog profile, which adjoins 
the Fermi rectangle at its upper edge. The total number of hard spheres in 
the Fermi rectangle should be the same as that of the added particles. One 
may obtain the same picture in a reverse way. Suppose we start from a high 
temperature where all the particles are active. We then slowly decrease the 
temperature to suppress the thermal motion. There will be a temperature 
where the density at the bottom is the closed packed density. Now, let us 
lower the temperature further. What happens? The next layer will become 
closed packed and is thus effectively frozen. As we keep lowering the tem- 
perature, the freezing of the particles will occur from the bottom and the 
frozen region will then spread out, until at T=0 all the particles are frozen. 




(15) 
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Note that the frozen particles in the closed packed region behave like a solid 
with <j) — 1. Such an observation helps us to resolve the puzzle associated 
with the disappearance of particles. The missing particles should form the 
condensate in the Fermi rectangle spanning the bottom to the lower part of 
the fiuidized layer. We term the surface that separates the frozen or a closed 
packed region with 0=1 from the fiuidized region with (j) < 1 the Fermi 
surface. The location of the Fermi surface, Cf, is determined by the number 
of the missing particles, namely, £p = p, — p / (3, and is thus a function of 
temperature. For nonzero we must put the missing particles below the 
Fermi surface and shift the bottom layer from £ = to From Eq.(15), 
we find that in 2D the number of missing or condensed particles, N a (T), at 
T < T c is N (T) = Q x (p - Ho/0) = N(l - T/T c ), where N(= Q x p) is the 
total number of particles, and T c is the condensation temperature defined as 
the point where the particle conservation, Eq.(15), breaks down, namely 

T c = mgDfi/fio (16) 
Hence, the fraction of condensed particles is given by: 

N (T)/N = 1 - T/T c (17) 

How is this picture modified in the presence of dissipation? With dissipation, 
solving the Enskog equation becomes a nontrivial task for three reasons: 
First, we do not yet know the precise functional form of the velocity distri- 
bution function for inelastic particles. Experimentally observed profiles [p3| 
definitely indicate that they are not Gaussian. Second, even if we somehow 
have empirical formulas for the velocity distribution functions, carrying out 
the integral with non-Gaussian profiles for the Enskog integro differential op- 
erator becomes a non-trivial task. Third, in the presence of dissipation, the 
temperature profile is nonuniform pij , and thus one has to solve the equation 



for the density profile, i.e., Eq.(8), along with the energy equation 25 . Ana- 
lytic solution for this case is difficult to obtain. However, if the dissipation is 
small, then one might assume that the velocity distribution function is still 
Gaussian and the temperature remains uniform except near the heat reser- 
voir. Under such assumption, we only need to solve the Eqs.(8)-(10) with the 
corrected pressure term Q, which is given by: 

P = pT[l + 1P D d X ] 

where 7 = ^(1 + e) when d — 2 and 7 = -|(1 + e) when d = 3. Since the 
solution of the force balance equation, dP/dz = —pg yields the same result, 
in the presence of dissipation, we use the force balance equation instead of 
solving the Enskog equation as was done in the elastic case eDefining 
cf> = pD d as above and using the above forms of P in the force balance 
equation yields the differential equation 

/30(C) + #(C)K = - 7 [^xW/# + 2 x (0)d0(C)/dC] (18) 
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Upon integration we find in the two dimensional case: 



P(C ~ P) = ~ log<£ + ci^ + c 2 log(l - a<f>) + c 3 /(l - a<f>) + c 4 /(l - a</>) 2 (19) 
where 

ci = -27a 2 /a 2 , 
c 2 = 7(ai/a 2 - 2a 2 /a 3 ) 

C3 = -C2 

c 4 = 7(— 1/a + ai/a 2 — a 2 /a 3 ) 

Then /3/i is the negative of the right hand side of (19) with cf) replaced with 
(j)o = 4>(( = 0). Integrating /3( from zero to (f>o yields 

(3/j, = c 5 + c 6 o - c 2 <t>l + c 7 / (1 - acf) ) + c 8 / (1 - acf) ) 2 (20) 

where 

c 5 = 7(— 2aia + a 2 + 3a 2 )/a 4 

C6 = -c 2 + 1 
C7 = — 7(2a 2 — 3aia + 4a 2 )/a 4 
c 8 = —a/ a 

This result may be shown to be equivalent to the elastic case when e = 1. 
Substituting the numeric values for the constants a, a\, a 2 and evaluating 
at 0o = 1 yields 

Aio = 1 + 10.3779(e+ 1) (21) 
In three dimensions we find 

Pit -fi) = - log* - - 1^^, (22) 

where f3ji is again given by the negative of the right hand side of (22) with 
*o replacing 4>. Then we integrate fi( between zero and 0o to get 

(Jfi — (po + -— -75 + 77: —To K Z6 > 

2(l-a0 o ) 2 2{l-acj> y 

where a = 77/6. When e = 1, this again reduces to the profile for the elastic 
hard spheres (12). We find upon evaluating (3/j, at <po = 1 that 

^ = l + 7.14964(e+l) (24) 

We note that in general, /j,q takes the form 1 + C(e +1), and we recall 
that T c = mgDfi/fio- Hence we may write 
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T c (e) = T c (e = 1) + - + — - ^ - ^ (25) 
Letting 5 = (1 — e)<cl, we may recast the above result as 



2 



T c (e) = T c ( £ = 1) + ( TTW ^ + ( TTW ^ + - (26) 

In summary, we conclude that the condensation phenomenon persists in the 
presence weak of dissipation. Note that in the above analysis, we set <f> = 1 
as the upper bound for a cage model. For maximum packing, <j) can exceed 
unity, as discussed before, i.(in 3D, <f> w 0.74 7r/6 and O = 2/y/3.) The clus- 
tering of particles near the bottom wall in two dimensional experiments p7[ 
may be a strong confirmation of this scenario. 



7 Conclusion 

In this review, we have explored simple consequences of excluded volume 
interaction for dense granular materials, and shown that the granular statis- 
tics cross over from Boltzmann to Fermi statistics as the strength of the 
excitation is reduced. We have also advanced a theory that such a crossover 
may be understood as the condensation of hard spheres below a condensa- 
tion temperature. Since the Enskog equation breaks down at high densities, 
comparable to the closed packed density, we have constructed physical solu- 
tions beyond this regime by forming the Fermi rectangle near the bottom. 
Those particles in the Fermi rectangle are assumed to be condensed. This 
way, we have extended the applicability of the Enskoq theory to arbitrarily 
low temperatures. 

As demonstrated in Eq.(14), the Enskog solution based on a virial expan- 
sion from the dilute to the dense case conserves particles only when we allow 
the density profile greater than the maximum closed packed density at the 
bottom, which is clearly unphysical. This is due to the fact that the equation 
of state, i.e. the pressure as function of density, is not known for all densities. 
If one would use the correct equation for the pressure (which is not ana- 
lytically known at this point, except for a numerically determined one p9[) 
that must diverge at the close packed density, it may be possible that the 
strange condensation phenomenon associated with the Enskog equation may 
be absorbed by an advanced theory |^{|. We point out, however, that even 
with this empirical pressure profile the formation of Fermi rectangle is clearly 
visible. [p9| The density profile above the Fermi rectangle was shown to be 
nicely fit by the Enskoq profile. Thus, our picture of condensation, which 
is an alternative way to extend Enskog theory results up to the maximum 
possible density, seems to capture the essence of physics - and we have shown 
analogies to the Fermi statistics. 
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For a polydisperse system, the upper limit for 77 is usually larger than 
that of the monodisperse system. Thus, we expect the freezing temperature 
to be much lower for a polydisperse system. Note, however, that as long as 
the upper limit for r\ is smaller than unity, Enskog theory for a gas predicts 
the condensation phenomenon, in the sense that unphysical densities are 
predicted. For the extreme limit of an Appolonian packing p8| , for which the 
close packed density can be unity, 77 = a<f) in (14) can be arbitrarily close 
to unity and thus the right hand side of (14) can be arbitrarily large. In this 
case, the condensation phenomenon discussed in this paper disappears. 

The original experiment done in Ref. || used a rectangular box with tilted 
side walls, but one may use a different shape such as a parabola to mimic the 
density profile of the electron gas in 3D, which may help one to visualize the 
excitation spectrum near the Fermi surface. Finally, we point out that while 
it is trivial to drive hard sphere grains by a thermal reservoir in MD simula- 
tions, it becomes a nontrivial task to do this experimentally. The normal way 
of exciting the granular system experimentally is through vibration, yet it 
remains an unresolved issue to determine the precise relation between the vi- 
bration strength and the effective temperature of the reservoir. Eq.(5), which 
is based on a single ball picture, may be a first step in this direction. The 
discovery of the Fermi statistics and the associated condensation phenomena 
may provide an avenue to apply the methods developed in equilibrium statis- 
tical mechanics to the study of granular dynamics, notably from the point of 
view of elementary excitations such as the Fermi liquid theory in condensed 
matter physics. 
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Density Profile for 100 Particles, f}=4, T=2.0 and T=20.0 
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